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A coupled map model for the chaotic phase synchronization and its desynchronization 
phenomenon is proposed. The model is constructed by integrating the coupled kicked oscil- 
lator system, kicking strength depending on the complex state variables. It is shown that the 
proposed model clearly exhibits the chaotic phase synchronization phenomenon. Further- 
more, we numerically prove that in the region where the phase synchronization is weakly 
broken, the anomalous scaling of the phase difference rotation number is observed. This 
proves that the present model belongs to the same universality class found by Pikovsky et 
al.. Furthermore, the phase diffusion coefficient in the de-synchronization state is analyzed. 
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q \ §1. Introduction 



Coupled chaos systems are not only theoretically interesting but also important 
from the viewpoint of the application to engineering and biology. Dynamical 
behaviors observed in coupled oscillator systems depend on the number of chaotic 
oscillators as well as the coupling form. As the number of oscillators is increased, 
the variety and the complexity of dynamical behaviors increase. One of the eminent 
characteristics of the coupled chaos systems is the synchronization and its break- 
down. For a coupled system consisted of identical chaotic oscillators, it was shown 
that the oscillators can completely synchronize under certain conditions. 4 ) This phe- 
^ ■ nomenon is called either the chaotic complete synchronization or the chaotic Huygens 

phenomenon. 4 ) On the other hand, even when two oscillators have a mismatch in 
characteristic, they can show the so called chaotic phase synchronization for a certain 
region of coupling strength. 5 ) 

It is well known that the coupled map system played a significant role in the 
study of complete synchronization. 7 )' 8 ) Particularly, coupled one-dimensional map 
system is a typical model in discussing the chaos synchronization-desynchronization 
phenomenon. In contrast to the complete synchronization, coupled mapping model 
for chaotic phase synchronization-desynchronization phenomenon is not studied so 



*' Email: fujisakaOi. kyoto-u.ac.jp 
**' Email: uchiyama@amath.hiroshima-u.ac.jp 

***) Present address: Department of Mathematical Sciences, Osaka Prefecture University, Sakai, 
Osaka 599-8531, Japan, Email: horita@ms.osakafu-u.ac.jp 



typeset using VTpT^.ds (Ver.0.9) 



2 



H. Fujisaka, S. Uchiyama, and T. Horita 



much because of the lack of suitable mapping model for coupled phase synchroniza- 
tion. The aims of the present paper are to propose a coupled mapping model for 
chaotic phase synchronization and then to examine the universality of the breakdown 
of synchronization so far mainly studied for the coupled Flossier system. 

The present paper is constructed as follows. In Sec. 2, we propose a general, 
linearly coupled mapping system for coupled chaotic oscillators which show either 
periodic or chaotic characteristics. The model is constructed by integrating the 
equations of motion for oscillator variables with a state dependent kicking term. In 
Sec. 3, using a two oscillators system, we discuss general characteristics of the system, 
and then study the phase difference statistics associated with the breakdown of the 
chaotic phase synchronization. We give summary and concluding remarks in Sec. 4. 

§2. Mapping model of coupled oscillatory chaotic systems 

First, consider the linear equations of motion 

m= P , (2-i) 

rp(t) = icuip - p, (2-2) 

where ip and p are complex variables, r is a positive number and to is a real number. 
In the limit r — > 0, the above equations reduce to ip = iujip, the equation of motion for 
a harmonic oscillator. So, u has the meaning of the eigenfrequency of the oscillator, 
and the above set of equations of motion turns out to describe a harmonic oscillator 
slightly modulated by the introduction of the "inertia term" rp. The above equations 
of motion can be solved as ip(t) ~ e xt , p(t) ~ e xt , with the characteristic exponents 
A = t -1 (— 1 ± y/1 + A%tuj)/2. For small r, they are written as A = iuj + tu 2 , — r _1 — 
ioj. One should note that in the presence of the inertia term, the above harmonic 
oscillation is unstable. 

We consider the periodically-kicked harmonic oscillator system 

m= P , (2-3) 

+oo 

rp{t) = i^- P +F a {^,r)- <*(*-*»)> (2-4) 

n=— oo 

(t n = nT), where T is the period of kicking and is taken to be unity without loss 
of generality. F a (ip,ip*) is the state-dependent complex amplitude of kicking, a 
standing for a set of parameters which characterizes the oscillator. Let us solve the 
above equation for t n — 5 < t < t n+ \ — 5, 5 being a positive infinitesimal quantity. 
Then, taking the limit r — ► 0, one finds that the above set of equations of motion 
reduces to 

^ n+1 = e^f a (^ n ,r n ) (2-5) 

and p(t n +i — 5) = iujip n+ i, where we noticed that ip n = ip(t n — 5) = ip(t n ) because 
of the continuity of if>(t) at t n . Here we defined 



(2-6) 
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Equation (2 • 5) with eq. (2 • 6) is the mapping system for the kicked-oscillator system 
(2 • 3) and (2-4). One should note that the mapping model (2 • 5) has two types 
of characteristic parameters. One is the eigenfrequency u> of the oscillation and 
therefore is relevant to the phase degree of freedom, and the other is the parameter 
set a, which is regarded as the control parameter mainly relevant to the amplitude 
dynamics of oscillation. 

In the above, we proposed a mapping model relevant to a chaotic dynamics show- 
ing a well-defined oscillation. Next, we will construct a mapping system consisted of 
coupled oscillators. As an example, consider the two oscillators system coupled to 
each other, 



7 pU) =p U) i (2-7) 
rp® = iurfU) - p U) + F aj {^\ V 0> ) • f t n ) + — - (2-8) 



for (j, k) = (1, 2), (2, 1). The last term of eq. (2 • 8) represents the coupling and K 
is the coupling constant. In the present paper, unless it is stated, K is kept to be 
non-negative. For the coupled iV-elements system, the equations of motion for the 
j-th oscillator are written as 



^=p^\ (2-9) 

+oo 

Tp® = iuj^ {j) - p (j) + F aj (V> (i) , ^ 0> ) • ^2 S(t- tn) + V^ j \ (2-10) 



Here T>tp^ represents the linear coupling term. Particularly, for the N = 2 model 
given in eqs. (2 • 7) and (2 • 8), we get Z>V (1 ' 2) = (^/2)0 (2,1) - V> (1 ' 2) )- 

Solving the above equations of motion for t n — 5 < t < t n+ \ — 5 and taking the 
limit r — > 0, one obtains the coupled map system 

= J2 J ^ U ' k faM ) ^ { n > ), (2-H) 
k 

for = ^\t n -5) = ^\t n ) <mdp^(t n+l -5) = (ico, +V)^ ] +1 , where / (^*) 
is the same as in eq. (2-6). The coupling constant Jjk has been defined as follows. 
Let us introduce the matrix A with the jk element Ajk defined by 

{iuj j + V)g j = J2 A jk9k (2-12) 

k 

for an arbitrary quantity Qj defined for the oscillator j. The interaction kernel Jjk 
in eq. (2-11) is defined via 

e iu j+ v g . = J2 (V) . k 9k = Y. J 3k^ k 9k- (2-13) 
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Equation (2 • 11) together with eqs. (2 • 12) and (2 • 13) is the fundamental result of 
the present paper. 

Particularly, if the coupling term has the form 



N 



Vgj = Yj K Jk(.9k -9j), 



(2-14) 



k=l 



where K~ k is the coupling constant, then the matrix element is given by 



A jk = iujj5 jk + K jk - K jij S jk- 



(2-15) 



If all characteristic frequencies are same, i.e., cj± = o>2 = ■ • ■ = con = w, and the 
coupling operator satisfies the condition ^ ■ T>gj = and therefore ^ ■ Jj k = 1, then 
eq. (2 • 11) is rewritten as 



(J) 

n+1 



^ faM^^*) + E J ^ [-U^V^*) " /o,^,^*)! (2-16) 



Furthermore, for the global coupling Vgj = (K/N) J2k=i(dk ~ 9j), one obtains Jj k = 
K/N for any combinations of j and k. In addition, when the system characteristics 
are all the same (a\ = a 2 = ■ ■ ■ = a,N = a), one finds that the equations of motion 
under consideration has the complete synchronization tpn^ = = • • • = = 
V>°, which obeys V'n+i = e *^/a(V'n' V'n*) irrespectively of the coupling strength. This 
fact suggests the possibility of the existence of the transition between the complete 
synchronization and its broken state. Details will be given in Sec. 3 for the two- 
oscillators system. 

§3. Chaotic phase synchronization in a two oscillators system 



In this section, we study the synchronization and desynchronization for the cou- 
pled maps model of a two oscillators system proposed in the preceding section. For 
the coupling term Vg lj2 = (K/2)(g 2 ,i - 51,2), we get 



A 



iuj\ — y 

K ■ ~ K 

T 2 ~ T 



K 
2 



(3-1) 



This immediately gives the coupled map system 



m - J K (u! - W2 )e^/a 1 (^ 1) ,4 1) *) + J' K ("i -U2)e^f a2 {^\^>), (3-2) 

Chi = J'k(W - ^e^faMl^l + Jk(u 2 - U^fa^^l, (3"3) 



V'n+l 
,(2) 



where the coupling constants are 
J k (Auj 



, V « 2 - (^) 2 A smh^S 
cosh — + iAuj- 



^K 2 - (Alo) 2 



(3-4) 
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, , * ^ i a, k K smh tt — — 

It should be noted that the coupling constants satisfy Jk(Alo)* = Jr{— Auj) and 
J' K (Au)* = J'A-Au). 

3.1. General characteristics of the two oscillators system 

Since Jq(Auj) = 1 and J' G (Au) = 0, oscillators are independent of each other for 
K = 0. On the other hand, for K — >■ oo, noting J^Aj) = ie"^^ = ./^(Aj)*, we 
obtain 

V4 1} = V4 2) = (3-6) 
= ±e*<"*+<«> [f a Ml€*) + fa 2 {€,€*)] • (3-7) 

This result implies that even if the characteristics of two oscillators are different, 
the system shows a complete synchronization in the strong coupling limit K — ► oo. 
It should be noted that for the coupled system consisted of identical oscillators the 
complete synchronization is achieved irrespectively of the strength as a particular 
motion, while the present synchronization for non-identical oscillators is observed 
only when the K — ► oo limit. *) On the other hand, if the system parameters Uj 
and aj are different from each other, the system does not show a complete synchro- 
nization. Nevertheless, we expect that there exists a certain kind of synchronization 
state if the coupling constant K is suitably chosen, which is nothing but the phase 
synchronization. 5 ^ This will be studied in the following subsection. 

Particularly, if u\ = uji = to and a\ = ai = a, i. e. , the two oscillators are 
identical, then the coupled map system is written as 

rfi = e*" + 4 [/a(Y4 2) ,V4 2) *) - /a(l#W)] } , (3-8) 

Y& = {U^\^l + 4 [faW>P,W) - faW, Vi 2) *)] } , (3-9) 

where J' K = (1 - e~ x ) /2. In this case, the particular, complete synchronization state 
exists and the dynamical variable Y>°(= = V^) obeys Y'n+i = e luJ fa(ipni "0n*)- 
Around the complete synchronization state, the small deviations <5^4^ = — V'n 
and <5V4 2 ^ = i^i? — i^n from the complete synchronization obey 

= e iu! {G lin 6iP+ + G 2 , n #+*), (3-10) 
^"+1 = e^e- K (G hn 5ii- + G 2 ^~*), (341) 



*•* The above result holds for K > 0. In the opposite limit (K — > -co), we observe the "anti- 
phase" oscillation Vn } = -V>£ 2) . Namely for if -> -oo, noting J K (Z\w) = \e~i AoJ+m , J' k {Alo) = 
— Jk{Au>)* , one gets the equation of motion 
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where ty± = (^ , ±^ , )/2, and Gi, n = df a (^ n ,^ n *)/d^ n and G 2 , n = df a {^l*)/d^\ 
It should be noted that the two equations of motion for 5tp^ are separated from 
each other, and the first equation is identical to the perturbation equation <5^° +1 = 
e tuJ (Gi jn 8ijj^ l + G 2<n 5ijj®*) f° r the small change of the initial condition of the element 
dynamics. Let A be the largest Lyapunov exponent of the element dynamics. One 
finds that 

|^+|~e A "|#+| (3-12) 

for large n. The exponent A shows the characteristic of the uncoupled system and 
is free from the stability of the complete synchronization. On the other hand, the 
second equation gives 

\w\ = l\i>P-i>P\~^ ±n m\, (3-i3) 

where 



X± = X — K. (3-14) 

The parameter Aj_ called the transverse Lyapunov exponent or the stability parameter 
determines the stability of the complete synchronization. 4 ) Namely, if Aj_ < 0, the 
complete synchronization is linearly stable. On the other hand, if A^ > 0, the 
complete synchronization is unstable. If the element dynamics is periodic, then 
A < and therefore Aj_ < 0, which implies that the complete synchronization is 
linearly stable. On the other hand, for a chaotic element dynamics (A > 0), the 
system shows the transition at K = A(= K c ) as the coupling strength K is decreased 
from above K c , and one observes the on-off (modulated) intermittency for K slightly 
below if c . 9 )< 10 ) 

3.2. Chaotic phase synchronization and its breakdown 

In the present paper, as a model of element dynamics, we use the mapping 
function 

f a (iP,r) = [a-(l + ibM 2 ]iP, (3-15) 

(i) 

where the parameters a and b are chosen in such a way that the dynamics V'n+i = 

z 1 ^ 3 faj{i>n\i>rP*) shows a chaotic behavior. One should note that the statistical 
characteristic of the isolated element dynamics is free from the choice of Uj because 
the local dynamics does not depend on the eigenfrequency after the replacement 
ipn^ —> e luJjn ipn ^ ■ A trajectory of the element dynamics (3- 15) for to = 0.08, a = 2.55 
and b = is shown in Fig. 1. Figure 2 depicts how the trajectories depend on b. One 
observes that depending on b, the isolated dynamics shows various characteristics. 
In later numerical experiments of coupled map system, we use the same values for 
a and b as a\ = a 2 = 2.55 and b\ = b 2 = for which the element dynamics are 
chaotic, where the largest Lyapunov exponent A is about 0.144, and the frequency 
parameters uo\ and u> 2 are chosen to be set as Auj = u\ — uj 2 = 0.08 except the figure 
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4. One should note that in the present model only the difference of eigenfrequencies 
is relevant. It was numerically proved that as far as b is small enough, the intro- 
duction of non-vanishing bj does not change the qualitative results found for bj = 0. 



2 




-2-10 1 2 



ReM 

Fig. 1. Trajectory for the element dynamics 
(2 • 5) with the mapping function (3 • 15), 
where w = 0.03, a = 2.55 and b = 0. The 
dynamics is chaotic with the largest Lya- 
punov exponent A = 0.144. 



Let 9n be the phase of ip„ r , i.e., 

^ = I^V^- It should be noted 
that as far as the phase variation is small 
enough, the phase dynamics of 9){ is 
uniquely determined even for a coupled 
system. For the control parameters aj 
and bj used in the present paper, the 
phase variations are small enough. The 
phase difference A9 n = 0$ — 9n^ obeys 
the equation of motion 

A9 n+1 = A9 n + a n (A9 n ), (3-16) 

where a n (A9), which is uniquely deter- 
mined by eqs. (3 • 2) and (3 • 3), is a 2tt- 
periodic function of A9, i.e., (a n (A9 + 
2ir) = a n (A9)). The integration of the 
above equation yields 



n-1 

A9 n = A9 Q +Y, a m {A9 m ). (3-17) 

m=0 

By taking the statistical average of the above equation, which is equivalent to the 
time average with the ergodicity assumption, the angular velocity of the phase dif- 
ference (rotation number) is given by 

1 n— 1 

AQ = hm - V a m (A9 m ) = (a), (3-18) 

n— *oo 77, * — » 

m=0 

where (• • • ) is the statistical average over an ensemble representing the steady state 
of the present chaos. Therefore, the average value of the phase difference obeys 

(A9 n ) = A9 + AQ ■ n. (3-19) 

Figure 3 shows temporal evolutions of the phase difference for different values of 
the coupling strength K. One observes that depending on K, the coupled mapping 
system exhibits chaotic phase synchronization {AQ = 0). Figure 4 displays how 
the rotation number AQ depends on K and Au(= u\ — W2). One finds that under 
the change of either K or Aoo a clear transition between phase synchronization and 
phase desynchronization state occurs. 

In order to study the critical behavior of AQ as the coupling strength K is 
decreased by keeping Aoj constant to be 0.08, the quantity (AQ) 2 is plotted as a 
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Fig. 2. Trajectories for the element dynamics (2 • 5) with the mapping function (3 • 15) for a = 
2.55 and various values of b. The abscissa and the ordinate are respectively Re(e" m ^i„) and 
Im(e _! ^t/; n ). Trajectories do not depend on the value of lo for the mapping (3 ■ 15). 



function of K in Fig. 5(a). The figure suggests that the critical behavior 

AQ ~ (K* - K)* 



(3-20) 



holds in the region K < K*, where AT*(« 0.0612) is a characteristic value of K being 
below K c (m 0.0804), the desynchronization point, i.e., AQ = for K > K c . The 
critical behavior (3 • 20) is identical to that of the saddle-node bifurcation, and is 
well-known in association with the existence of the type I intermittency. Figure 5(b) 
shows the numerical result for the K region below the onset point K c of the chaotic 
phase synchronization. Slightly below K c , the rotation number turns out to take the 
asymptotic form 



AQ ~ exp -c • (K c - K)~ 



(3-21) 



instead of the saddle-node type, where c is a positive constant. This behavior is a 
result of chaotic, one may say, "stochastic" emergence of the saddle-node type chan- 
nel, and has been reported near the desynchronization point in a phenomenological 
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Fig. 3. Temporal evolution of phase differ- 
ences for the coupled two maps system. For 
K = 0.081, the phase synchronization is 
achieved. 



Fig. 4. The rotation number AQ depends on 
both the eigenfrequency mismatch Au)(= 
oji — LU2) and the coupling constant K. 
For AQ — 0, the phase synchronization is 
achieved. 



mapping model 13 -* and experimentally in the modulated CO2 laser system. 12 ) So, the 
present model belongs to the same universality class of chaotic phase synchronization 
showing the above two different critical behaviors. 
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K (K c -K)- 1 ' 2 
Fig. 5. Critical behaviors of the rotation number Af2 as functions of the coupling constant K. 

The fluctuation of the phase difference A6 n from the average value (3 • 19) is 
evaluated with the variance 

a 2 n ^{{Ae n -{Ae n )f). (3-22) 

It is expected that by assuming the mixing property of a n (A8 n ), the variance obeys 
the diffusion law 

ct 2 = 2Dn (3-23) 

for large n, where the diffusion coefficient is given by D = (1/2) X)nL-oo CLi with 
C n = {5a n (A6 n )5ao(A9o)) , (5a n = a n — (a)). In the present paper, the diffusion 
coefficients were calculated with the asymptotic form (3.23) for various values of 
K by keeping Aoo = 0.08. Results are shown in Fig. 6. One observes an anoma- 
lous enhancement of the diffusion coefficient near K*, the lower synchronization- 
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desynchronization transition point. Figure 7 displays the relation between the rota- 
tion number and the diffusion coefficient. One observes a linear relation D ~ AQ for 
K slightly below K c . This fact agrees with the observation in the coupled Rossler 
system. 14 ) The numerical result shows D j AQ « 2.7. 




0.02 0.04 0.06 0.08 



K 

Fig. 6. The rotation number AQ and the 
phase diffusion coefficient D below the 
breakdown point K c of chaotic phase syn- 
chronization. The diffusion coefficient 
shows an enhancement slightly below the 
breakdown point of the phase synchroniza- 
tion. 



0.012 



0.004 
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An 

Fig. 7. The relation between the phase dif- 
fusion coefficient and the rotation number 
near the breakdown point of the chaotic 
phase synchronization. One observes a lin- 
ear relation D ~ AQ. 



§4. Summary and concluding remarks 

In the present paper, we proposed a mapping model of coupled chaotic oscilla- 
tor system composed of element dynamics each of which has a periodic oscillation 
characteristic in an uncoupled case. An oscillation has two degrees of freedom, phase 
and amplitude. Parameters contained in dynamical systems described by differential 
equations of motion for describing oscillatory behaviors generically control both the 
phase and the amplitude simultaneously. On the other hand, the present mapping 
model, the parameter u mainly controls the phase dynamics and the parameter a the 
amplitude dynamics. Therefore, in order to study the characteristic dynamics of the 
phase and the amplitude separately, the present mapping model of oscillatory dy- 
namics is considered to be more convenient to analyze chaotic phase synchronization 
than the differential equation system such as the Rossler model. 

By making use of concrete coupled two maps system, it was shown that the 
present coupled map system shows the phase synchronization-desynchronization phe- 
nomenon. We found that when the amplitude synchronization is achieved, after the 
breakdown of the phase synchronization the rotation number turns out to show a 
normal scaling (3 • 20). However, since this characteristic is quite particular to the 
element map used in the present paper, in general no amplitude synchronization is 
expected to be observed, and it may be concluded that sufficiently near the tran- 
sition point the anomalous scaling (3 • 21) is ubiquitously observed. In this sense 
the present two maps system belongs to the universality class of that reported in 
Ref. 13). Furthermore, we studied the phase diffusion coefficient, and observed that 
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its critical behaviors are the same as those found in Ref. 14). 

So far, mapping models of coupled chaos systems, particularly coupled one- 
dimensional mapping models, have been extensively studied to clarify their varieties 
of dynamical phenomena as well as mathematical structures of their dynamical be- 
haviors. The proposed coupled map model composed of chaotic oscillator elements 
which show well-defined limit-cycle type oscillation characteristics is expected to con- 
tribute to the progress in the study of the phase synchronization-desynchronization 
phenomenon in a wide class of coupled chaos systems. 
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